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Abstract. In this study we investigate how the condition of zero 
current on open flux tubes with polar wind outflow, subjected to 
large photoelectron fluxes, can be achieved. We employ a steady 
state collisionless semikinetic model to determine the density 
profiles of 0 + , H + , thermal electrons and photoelectrons coming 
from the ionosphere along with H + , ions and electrons coming 
from the magnetosphere. The model solution attains a potential 
distribution which both satisfies the condition of charge neutrality 
and zero current. For the range of parameters considered in this 
study we find that a 45-60 volt discontinuous potential drop may 
develop to reflect most of the photoelectrons back toward the 
ionosphere. This develops because the downward flux of 
electrons from the magnetosphere to the ionosphere on typical 
open flux tubes (e.g. the polar rain) appears to be insufficient to 
balance the photoelectron flux from the ionosphere. 


Introduction 

Photoelectrons are produced in the upper atmosphere during 
the ionization of neutral species by solar EUV radiation. They 
have long been suspected of playing a role in determining or 
moderating the polar wind [ Axford , 1968; Lemaire and Scherer , 
1972]. Because of their high energies (~ 20 eV), many of them 
are able to freely flow along magnetic field lines out of the 
ionosphere. The typical escape flux of these electrons ranges 
from 10 9 cm' 2 s' 1 (at solar minimum) to 6 x 10 9 cm' 2 s' 1 (at solar 
maximum). These values come from a two stream phototelectron 
model [Richards et al , 1994] which includes all of the relevant 
ionospheric processes involved in the production and transport of 
the photoelectrons, and has been extensively calibrated against 
AE and DE photoelectron measurements. Such fluxes represent a 
current density, at the source altitude, of 1 .6 pA/m 2 to 9.6 pA/m 2 
respectively. These currents are comparable to those seen on 
field lines in the auroral zone. Since many noncurrent-carrying 
field lines have such photoelectron fluxes on them the plasma 
must adjust itself in some way so that a large current does not 
develop. How is zero current flow achieved? 

On closed flux tubes zero current could be achieved by a 
balance between outflowing photoelectrons and photoelectrons 
from the conjugate ionosphere. If such conjugate fluxes are 
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unequal then thermal electrons might be forced to flow in the 
opposite direction from one ionosphere to the other. Such forced 
inflow of the thermal electrons could heat them, due to adiabatic 
compression, and modify the ionosphere to generate ion 
outflows, if compressional heating dominated heatflow. 

On open field lines, how the plasma distribution achieves zero 
current is a more difficult question. Several possibilities may be 
proposed: 

1. The ions in the ionosphere could be accelerated outward by 
a potential drop so that they produce upflowing fluxes 
comparable to those of the photoelectrons. Since these fluxes 
exceed the hydrogen limiting flux the 0 + ions would be required 
to carry the bulk of such ion flux. This appears to be the 
conceptual nature of the solution in Tam et al. [1995]. 

2. High altitude thermal electrons (< few eV) could be drawn 
down the field line as in the hypothesized closed field line case. 
It is unknown if such electrons exits in large numbers. 

3. Magnetospheric electrons (> few eV) could be drawn down 
the field line to balance the photoelectron flux. Typical polar rain 
fluxes in the topside ionosphere are only 10 7 - 10 8 cm' 2 s' 1 
however, ( Winningham and Heikkila [1974]) which appears to be 
insufficient. 

4. A potential barrier could develop that reflects most of the 
photoelectrons back to the ionosphere. This, in principle, could 
be in the form of a narrow shock like structure [Barakat and 
Sc hunk, 1984]. 

Evidence Supporting the Existence of a Potential 
Barrier 

We explore here the last of these possibilities, namely that 
zero current is primarily achieved on open polar cap field lines 
through the development of a potential barrier that prevents the 
outflow of most of the photoelectrons. This is supported by 
LAPI observations of downflowing photoelectrons with energies 
up to 60 eV (or more) seen in the polar cap at altitudes below 
1000 km ( Winningham and Gurgiolo [1982]; Pollock et al. 
[1991]; Horwitz et al. [1992]), and by the observations of high 
speed ion beams seen on dayside polar cap field lines at high 
altitudes by the RIMS instrument on DE-1 ( Pollock et al. [1991]; 
Horwitz et al. [1992]). The stronger of these two pieces of 
evidence is the observation of reflected photoelectrons since the 
relationship between the ion beams and the photoelectrons is not 
conclusively established. 

If such a potential barrier forms and may be considered 
localized, where would it be expected to develop? This barrier 
likely forms at altitudes well above the topside ionosphere for the 
following reasons. 1) The average flow speeds of the 0 + and H + 
ions in the polar cap at altitudes <1 R E do not show any evidence 
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of large potential differences (5 volts or more) between the F 
region and an altitude of 1 R E {Chandler et al [1991]; Abe et al 
[1993]). 2) A collisional, semikinetic, transition region model 
(Wilson [1992], Wilson et al. [1996]), which does not include any 
flow enhancing effects that might be attributable to 
photoelectrons, predicts H + flow speeds that are about a factor of 
2 larger than the average observed below 1 R E altitude. Since 
results from this model are not driven by the boundary conditions 
and the model is not subject to sonic or heatflow singularities as 
are steady state generalized transport models, the difference 
between model results and data strongly suggest the operation of 
an additional drag force acting on the H + . If photoelectrons are a 
strong accelerator of outflowing ions the disparity between 
model results and data is greatly compounded. 3) The 
downflowing photoelectrons seen by the DE2/LAPI instrument 
show highly spatially structured characteristics even when the 
upflowing photoelectrons are stable over large distances 
[Winningham and Gurgiolo, 1982]. This may indicate the larger 
spatial structure of the distant magnetospheric region where these 
electrons are reflected or it may indicate the highly dynamic or 
unstable nature of the reflection process. 

Description of the Model 

The results presented in this paper were found with a steady 
state model which determines density profiles of the various 
species using the collisionless semikinetic treatment. In this 


method one can find an expression for the density of a given 
species as a function of position along the field line and the local 
value of the electric potential. The expressions derived for this 
paper assumed that the potential was zero at the lower boundary 
and decreased monotonically to the upper boundary (i.e. no 
trapping zones). The velocity distribution of all species at their 
respective boundaries for in going particles was assumed to be 
the appropriate half of an isotropic Maxwellian. The one 
exception was the photoelectron distribution which had a cutoff 
at 60 eV [Lee et al . , 1980]. The ionospheric thermal electrons 
were assume to be Boltzmann distributed. The assumption of 
Maxwellian velocity distributions for particles entering the model 
domain does not mean that the velocity distributions remained 
Maxwellian everywhere. The effects of magnetic folding, 
electric field acceleration and the loss cone combine to create 
fairly nonMaxwellian velocity distributions at points along the 
field line. 

The model assumes that all species are coliisionless. It can be 
shown from a mean free path analysis that 97% of the 
photoelectron flux is carried by photoelectrons with sufficient 
energy to be collisionsless given the densities and temperatures 
listed below. The thermal electrons will be heated by collisions 
with the photoelectrons but because of their high thermal 
conductivity the thermal electron temperature will remain within 
a few percent of being constant and isotropic. The 0 + ion flux 
and density profile will be unaffected by self collisions but these 
ions will be slightly heated by collisions with the thermal 





Figure 1 . Solar maximum case (with Dhotoelectron upward flux = 6 x 10 9 cm' 2 s’ 1 at 500 km) when the mantle plasma density at 9 R F 
is 0.05 cm-- 1 . The 0 + flux is 3.4 x 10® cm' 2 s' 1 , the H + flux is 2.9 x 10 7 cm' 2 s' 1 , the photoelectron flux is 5. 1 x 10 7 cm' 2 s'' 1 and the 
polar rain flux is 1.3 xlO 7 cm' 2 s' 1 . The combination of these fluxes leaves a residual current of 0.0089 pA/m 2 into the ionosphere. 
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electrons. This effect is accounted for through adjustment of the 
boundary 0 + temperature. The one species significantly affected 
by low altitude collisions is H + originating from the ionosphere. 
To find its density profile we have joined a low altitude 
collisional semikinetic model profile to that produced by the 
collisionless semikinetic model, at an altitude where collisions 
are infrequent. 

The parameters used in the cases shown in this paper were: 1 ) 
For ionospheric H + and 0 + , T 0 = 2000 K, N Q (H + ) = 1000 cm' 3 , 
and N 0 (0 + ) = 60000 cm' 3 . 2) For the thermal electrons T e = 
5000 K. 3) For magnetospheric electrons and ions, kT 0 = 100 eV 
with N 0 ranging between 0.5 and 0.05 cm' 3 . 4) For the 

photoelectrons, kT pe = 20 eV with N 0 ranging between 16.7 and 
100.2 cm -3 . These values of N 0 for the photoelectrons gave 
upward fluxes, in the absences of any electric potential, of 10 9 to 
6 xIO 9 cm' 2 s' 1 representing solar minimum to solar maximum 
conditions. 

To obtain a self consistent solution satisfying zero current and 
quasineutrality, the code iteratively modifies the electric potential 
until charge neutrality is achieved. The zero current condition is 
not directly applied to the iterative calculations but is only 
incorporated after a solution is found. Because of the 
transcendental nature of the model, multiple solutions are 
possible. One such solution involves a small potential drop ( < a 
few volts) along the flux tube, with the densities of the 
ionospheric ions and electrons equal, and separately the 
magnetospheric ion and electron densities equal. Such a solution 


does not produce a zero current because the potential drop along 
the field line is insufficient to reflect many photoelectrons back to 
the ionosphere. Another solution produces a large discontinuous 
drop in the potential. This feature is an electrostatic shock or 
rarefaction shock that can develop in systems with multiple 
electron populations that have significantly different temperatures 
(see Bezzerides et al. [1978] and Barakat and Schunk. [1984]). 
Such class of solutions can satisfy both the conditions of charge 
neutrality and zero current. 

Discussion 

Figures 1 and 2 show the results for two different cases with 
the photoelectron fluxes ranging between the solar minimum and 
solar maximum values and the magnetospheric plasma densities 
ranging between 0.05 and 0.5 cm' 3 . Such range of 
magnetospheric plasma densities are consistent with the plasma 
mantle [Shodhan et al . , 1996], but are low compared to the 
magnetosheath [Sibeck and Gosling , 1996]. They result in 
electron fluxes into the ionosphere of 10 7 and 10 8 cm' 2 s* 1 , 
consistent with polar rain measurements. In each figure the ion 
densities are plotted in the first panel, the electron densities in the 
second, the ionospheric H + and 0 + drift speeds in the third and 
the electric potential in the fourth panel. 

The total potential drop between the ionosphere and the 
magnetosphere in Figure 1 is 59.1 V compared to 47.4 V in 
Figure 2. This difference is due to the fact that the escaping 





Figure 2, Solar minimum case (with photoelectron upward flux = 10 9 cm' 2 s' 1 at 500 km) when the mantle plasma density at 9 R E is 
0.5 cm' 3 . The 0 + flux is 1.2 xIO 6 cm -2 s 1 , the H+ flux is 2.9 x 10 7 cm' 2 s' 1 , the photoelectron flux is 1.45 x 10 8 cm" 2 s' 1 and the polar 
rain flux is 1.23 x 10 8 cm' 2 s' 1 . The combination of these fluxes leaves a residual current of 0.013 flA/m 2 out of the ionosphere. 
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photoelectron flux is about 600 times greater than the 
precipitating magnetospheric electron flux for the case in Figure 
1 while it is only 10 times greater for the case in Figure 2. An 
expression that can be used to estimate the potential drop is given 
in equation (1) where the sum of the ion and electron fluxes is 
equated to zero (see Knight [1973]). 
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In this expression n eo and np eo are the magnetospheric and 
photoelectron densities at their sources, T tf and T are the 
magnetospheric and photoelectron temperatures respectively, 
A(j> is the potential drop (assumed to be positive), E m is the high 
energy cutoff for the photoelectrons (60 eV), and F, is the net ion 
flux along the flux tube. (In the absence of a solution one can use 
equation (1) to estimate the potential drop by assuming F i = 0.) 
This equation gives a potential drop of 46.9 V for the solar 
minimum case (Figure 2) and 59.2 V for the solar maximum case 
(Figure 1) using the ion fluxes found in those solutions. The 
results shown in Figures 1 and 2 leave a residual current density 
of -0.009 and 0.013 |iA/m 2 respectively at 500 km altitude. 

Unlike the results in Tam et al [ 1 995] (here after referred to as 
TM) the 0 + fluxes in these cases are relatively small; 3.4xl0 6 
cm' 2 s' 1 for Figure 1 and L2x 10 6 cm' 2 s’ 1 for Figure 2 at 500 km 
altitude. The oxygen ions do not participate much in zeroing the 
current. If the potential drop which allows attainment of near 
zero current were to develop at lower altitudes in the topside 
ionosphere it would not need to be as large as the results 
presented here since it would accelerated large fluxes of 0 + ions 
to escape energies so that these could contribute significantly to 
the current balance. This is the type of solution found in TM. 
There appears to be two primary choices for a solution that 
attains zero current on open, photoelectron carrying flux tubes. 
One is the TM result (and modified versions presented since) and 
the other is a solution with a high altitude strong photoelectron 
potential barrier, as presented here. Given that the TM result 
requires large, near limiting fluxes of 0 + to flow from the sunlit 
ionosphere and gives supersonic 0 + flows below 1000 km 
altitude (things which our solution does not give), we feel that 
our solution is closer to giving the actual effect of photoelectrons 
on geophysical polar cap flux tubes. 
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